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FOREWORD 


The  work  described  in  this  report  was  supported  by  the  U.  S. 

Air  Force  through  the  Hypersonic  Research  Laboratory,  Aerospace 
Research  Laboratories,  Office  of  Aerospace  Research,  Wright-Patterson 
AFB  Ohio,  under  Contract  F  33615-67-C- 1298.  The  contract  monitor  for 
this  program,  conducted  during  the  period  November  1966  through 
November  1969  was  Dr.  Wilbur  Hankey.  A  summary  of  the  studies 
conducted  under  this  contract  is  given  in  the  following  section.  A 
detailed  description  is  then  given  of  the  theoretical  and  experimental 
studies  of  shock  wave  -  laminar  boundary  layer  interaction  on  curved 
compression  surfaces  in  hypersonic  flow. 

The  author  wishes  to  acknowledge  the  contribution  of 
Mr.  J.  R.  Moselle  who  constructed  the  IBM  machine  code  used  in  the 
theoretical  study.  The  assistance  and  cooperation  of  Mr.  T.  J.  Bell  and 
Mr.  R.  F.  Hathaway,  of  the  Hypersonic  Facilities  Department  of  CAL, 
in  the  experimented,  program  is  also  gratefully  acknowledged. 
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ABSTRACT 


The  report  first  summarizes  the  studies  conducted  under  this  contract. 
A  detailed  description  is  then  given  of  the  theoretical  and  experimental  studies 
of  shock  wave-laminar  boundary  layer  interaction  on  curved  compression  sur¬ 
faces  in  hypersonic  flow.  In  the  theoretical  study  we  have  used  the  integral 
forms  of  the  boundary  layer  equations  for  the  conservation  of  mass,  stream- 
wise  momentum,  normal  momentum,  moment  of  streamwise  momentum  and 
energy  to  describe  the  development  of  attached  and  separated  boundary  layers 
on  a  curve  surface.  In  this  formulation  we  do  not  assume  the  pressure  to  be 
constant  across  the  boundary  layer.  The  theory,  which  was  in  good  agreement 
with  measurements  on  compression  surfaces  in  Mach  10  and  14.8  airflows, 
showed  that  large  pressure  gradients  can  be  generated  across  the  boundary 
layer  in  the  separation  and  reattachment  regions.  In  the  experimental  study, 
skin  friction,  heat  transfer  and  pressure  measurements  were  made  on  a  series 
of  flat  plate -cylindrical  arc-wedge  compression  surfaces.  Pitot  and  cone 
pressure  measurements  were  made  above  the  surface  in  the  separation  and 
reattachment  regions  of  separated  flows  to  estimate  the  static  pressure  dif¬ 
ference  across  the  boundary  layer.  The  radius  of  curvature  of  the  cylindrical 
arc,  the  inclination  of  the  wedge,  and  the  unit  Reynolds  number  of  the  free- 
stream  were  varied  to  examine  their  effect  on  the  properties  of  both  attached 
and  separated  interaction  regions  over  the  models.  We  found  that  separation 
first  occurred  on  the  compression  surface  downstream  of  the  flat  plate,  and 
incipient  separation  could  not  always  be  detected  by  observing  the  first  occur¬ 
rence  of  an  inflexion  point  in  the  pressure  distribution.  The  conditions  under 
which  incipient  separation  occurs  were  far  less  sensitive  to  surface  curvature 
than  those  observed  in  earlier  experimental  studies. 


A  SUMMARY  OF  THE  RESEARCH  PROGRAM 


The  research  program  which  has  been  completed  at  CAL  under 
Contract  No.  F  33615-67-C-1298  has  been  devoted  to  studies  of  the 
structure  and  development  of  attached  and  separated  boundary  layers  on 
compression  surfaces  in  hypersonic  flow.  The  emphasis  in  this  study  was 
placed  on  examining  shock  wave -laminar  boundary  layer  interaction  at 
free  stream  conditions  for  which  laminar  boundary  layers  could  be  ex¬ 
pected  to  be  encountered  on  high  performance  maneuverable  vehicles. 

Only  on  high  Mach  number  flight  (  Ma,  >  14)  at  high  altitudes  (low  Reynolds 
numbers)  on  highly  cooled  bodies  (tw/T0  ,  0  ( o.i )  )  are  laminar  boundary 
layers  of  practical  importance.  Under  these  circumstances  the  interaction 
between  the  viscous  and  inviscid  flow  over  the  leading  edge  (the  boundary 
layer  displacement  effect)  is  important;  and  since  the  leading  edge  must  be 
blunt  to  prevent  a  large  heating  load,  the  effect  of  the  highly  non-uniform 
inviscid  flow  field  (the  entropy  layer)  generated  by  bluntness  must  be 
examined.  In  many  cases  the  compression  surfaces  on  hypersonic  vehicles 
would  be  curved  and  flared  into  the  main  body  to  prevent  flow  separation, 
and  thus  we  set  out  to  find  how  surface  curvature  affects  the  development 
and  separation  of  the  boundary  layer,  and  how  large  are  the  normal 
pressure  gradients  developed  across  the  thick,  highly  curved  boundary 
layers  in  these  high  Mach  number  flows. 

A  paper  describing  the  theoretical  and  experimental  study  of 

leading  edge  bluntness  and  boundary  layer  displacement  effects  on 
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attached  and  separated  laminar  boundary  layers  was  presented  at  the 
AIAA  6th  Aerospace  Sciences  Meeting  in  New  York  (A1AA  Paper  68-68), 
and  is  to  be  published  in  two  parts  in  the  AIAA  journal.  A  paper  describing 
the  theoretical  and  experimental  studies  of  the  shock  wave -boundary  layer 
interaction  on  curved  compression  surfaces  was  presented  at  the  sympo¬ 
sium  on  "Viscous  Interaction  Phenomena  in  Supersonic  and  Hypersonic 
Flow"  sponsored  by  the  Hypersonic  Research  Laboratory  of  the  Aerospace 
Research  Laboratories,  Wright- Patterson  Air  Force  Base,  Ohio, 
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10-13  May  1969.  This  latter  work  is  described  in  full  in  this  report. 

Leading-Edge  Bluntness  and  Boundary  Layer  Displacement  Effects 

In  this  work  we  made  a  detailed  theoretical  and  experimental 
study  of  the  problem.  The  theoretical  analysis  used  the  integral  forms  of 
the  boundary  layer  equations  for  the  conservation  of  mass,  momentum, 
moment  of  momentum  and  energy  to  describe  the  development  of  the 
viscous  layer  and  its  interaction  with  the  outer  inviscid  flow.  The  analysis 
was  based  on  a  method  successfully  developed  earlier  to  describe  shock- 
induced  separated  flow  with  large  wall  cooling, but  applicable  only  to 
situations  where  the  interaction  upstream  of  the  corner  was  weak.  Here 
we  studied  flows  in  which  the  pressure  gradient  developed  by  the  viscous- 
inviscid  interaction  over  the  basic  flat  plate  was  strong  and  in  which  there 
is  a  strong  interaction  between  the  viscous  and  inviscid  flow  in  the  corner 
interaction  region. 

For  flows  over  highly  cooled  flat  plates  in  the  strong  interaction 
regime,  the  analysis  indicated  that  the  boundary  layer  became  supercriti¬ 
cal;  i.  e.  ,  it  will  thin  when  subjected  to  an  adverse  pressure  gradient.  In 
this  formulation  of  the  problem.in  order  to  join  the  solution  of  the  boundary 
layer  equation  close  to  the  separation  point  (where  the  boundary  layer  is 
subcritical)  to  the  integration  at  the  beginning  of  the  interaction,  it  is 
necessary  to  make  a  small  but  finite  jump  from  the  supercritical  to  the 
subcritical  branch  of  the  solution.  Here,  the  strength  of  this  jump  is 
specified  by  satisfying  the  conservation  of  mass,  momentum,  moment  of 
momentum  and  energy  across  the  jump.  Good  agreement  has  been  found 
between  theory  and  experiment  throughout  the  interaction  region  for  highly 

cooled  flows  in  the  weak  interaction  regime. 

« 

In  the  strong  interaction  regime  we  found  good  agreement 
between  theory  and  experiment.  However,  in  this  case,  we  experi¬ 
enced  difficulty  in  integrating  smoothly  through  the  neck  region  of  the  flow, 
and  we  had  to  introduce  a  scale  length  from  the  experimental  study  to 
specify  a  unique  solution.  At  this  time  it  was  believed  this  problem  arose 
because  the  pressure  gradient  developed  normal  to  the  viscous  layer  was 


v 


neglected  in  the  formulation  of  the  analysis.  Indeed,  this  was  found  to  be 
the  case,  for  in  the  analysis  which  we  developed  subsequently  in  which  the 
normal  pressure  gradient  was  taken  into  account,  we  experienced  no  pro¬ 
blems  in  integrating  smoothly  through  the  neck  region. 

The  experimental  investigation  of  le elding  edge  bluntness  and 
boundary  layer  displacement  effects  on  attached  and  separated  flows  over 
compression  surfaces  was  one  of  the  first  systematic  and  detailed  studies 
to  be  made  of  this  problem.  For  the  first  time  simultaneous  measurements 
of  skin  friction,  heat  transfer  and  pressure  were  made,  together  with 
schlieren  photographs  of  the  flows  on  a  series  of  sharp  and  blunt  flat 
plate-wedge  compression  surfaces  over  a  range  of  free  stream  Reynolds 
numbers  and  Mach  numbers.  As  in  earlier  studies,  we  investigated  the  time 
required  to  establish  a  steady  separated  flow  and  the  effect  of  a  finite 
model  span  on  the  size  and  properties  of  the  corner  interaction.  It  was 
found  that,  for  the  range  of  free -stream  conditions  and  model  configura¬ 
tions  studied,  attached  interaction  regions  were  established  within  the 
starting  time  of  the  tunnel  (less  than  1  millisecond).  However,  we  found 
that  separated  regions  could  take  as  long  as  2.  5  milliseconds  to  stabilize, 
and  the  mechanism  of  flow  establishment  was  consistent  with  the  propaga¬ 
tion  of  an  acoustic  disturbance  close  to  the  wall.  The  effect  of  finite  model 
span  on  the  properties  of  both  attached  and  separated  regions  was  investi¬ 
gated  by  varying  the  model  span  while  monitoring  both  the  longitudinal  and 
spanwise  properties  on  the  model.  Although  increasing  the  model  span 
from  1  ft  to  1. 5  ft  did  not  measurably  influence  the  properties  of  the 
attached  interaction  regions, it  did  cause  a  small  increase  in  the  size  of 
large  separated  regions.  However,  further  increasing  the  model  span  to 
2  ft  had  little  additional  effect,  and  the  addition  of  side  fences  to  this 
latter  configuration  did  not  affect  the  flow.  It  was  found  that  a  uniform 
distribution  of  properties  across  the  center  span  of  the  model  does  not 
necessarily  imply  that  the  measurements  are  free  from  span  effects. 

Also,  side  fences  should  be  used  with  care,  for  their  addition  to  some 
configurations  can  cause  as  much  flow  distortion  as  the  lack  of  them. 


In  the  experimental  studies  of  shock  wave -bound ary  layer  interaction 

on  the  highly  cooled  sharp  flat  plate-wedge  models,  which  was  performed 

at  Mach  numbers  from  14  to  20  and  unit  Reynolds  numbers  for  2200/in  to 

14,  000/in,  we  found  that  the  upstream  influence  on  the  wedge  in  attached 

flows  was  extremely  small.  However,  when  separation  is  promoted,  the 

interaction  feeds  forward  very  rapidly  as  the  wedge  angle  is  increased. 

The  abrupt  change  in  the  extent  of  the  upstream  influence  is  in  marked 

contrast  with  the  character  of  laminar  boundary  layers  in  supersonic  flow 

where  the  upstream  influence  increases  uniformly  with  wedge  angle  in  both 

attached  and  separated  flows.  The  length  of  the  separated  region  was  found 

to  decrease  with  increasing  free  stream  Mach  number  and  decreasing  unit 

Reynolds  number.  A  plateau  pressure  region  was  found  in  all  well 

separated  flows,  and  the  plateau  pressure  ratio  was  correlated  in  terms 

of  *  e  '/z  .  The  incipient  separation  condition  was  determined  at  each 

free  stream  condition  from  the  skin  friction,  heat  transfer  and  pressure 

measurements.  The  skin  friction  measurements  verified  the  separation 

criteria,  based  on  characteristic  changes  in  the  form  of  the  heat  transfer 

and  pressure  distributions  with  wedge  angle,  which  had  been  proposed  in 

earlier  studies.  We  found  that  those  measurements  of  the  incipient 

separation  condition  correlated  in  the  form  MS.  .  =4.  3*  '*  . 
r  mcip  *l 

The  study  of  the  interaction  regions  on  the  blunted  flat  plate -wedge 

was  performed  at  a  free  stream  Mach  number  of  19.  8  and  unit  Reynolds 

numbers  of  8700/in  and  2055/in  for  a  series  of  leading  edge  thicknesses 

(from  0.  007  to  1.  2  inches)  and  shapes  (cylindrical  and  flat)  to  obtain 

conditions  on  the  basic  flat  plate  from  the  viscous  to  the  bluntness- 

dominated  regime  (0.08  <  <  2).  The  measurements  of  skin  friction, 

heat  transfer,  and  pressure  made  in  the  attached  flow  over  the  blunted 
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flat  plates  agree  well  with  the  theory  of  Cheng  et  al  from  the  displacement  to 

the  bluntness-dominated  regime  (2  >  *  //f  >  0.  08).  For  well 

£  f  /Vi 

separated  flows,  increasing  the  leading  edge  bluntness  for  >  0.  5 

slightly  increased  the  size  of  the  separated  region,  whereas  for  Kf/Kj  c  0.  5, 
an  increase  in  the  bluntness  caused  a  marked  decrease  in  the  length  of  the 
separated  flow.  A  well  separated  region  in  these  flows  was  always 
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characterized  by  a  plateau  pressure  region.  The  plateau  pressure  ratio 
decreased  with  increased  bluntness,  and  it  was  found  that  the  measurements 
correlated  in  the  form  Pplateau /#,  (-^-f  versus  l(g/  K  .  The 

pressure,  heat  transfer,  and  skin  friction  on  the  wedge  in  both  attached 
and  separated  flows  are  markedly  reduced  by  leading  edge  bluntness. 

Thus,  in  the  design  of  compression  surfaces,  a  compromise  must  be 
reached  between  protecting  the  leading  edge  and  recovering  as  much 
pressure  as  possible  on  the  ramp.  We  found  that  in  flows  strongly 
influenced  by  leading  edge  bluntness,  the  skin  friction  measurements 
were  the  only  accurate  method  of  determining  incipient  separation,  though 
the  heat-transfer  measurements  still  provided  a  qualitative  separation 
criterion.  The  wedge  angle  required  to  cause  incipient  separation  was 
reduced  by  leading  edge  bluntness,  though  the  trend  was  weak;  the  measure¬ 
ments  correlated  in  the  form  M6  .  .  Apt-)  -  ¥-.3f  /-— *777 ) 

mcip^,L/  Uf  / 

A  Solution  to  the  Second  Order  Boundary  Layer  Equations  Including  Normal 
Pressure  Gradient 

Ir.  order  to  describe  properly  the  viscous -inviscid  interaction  in  the 
separation  and  reattachment  regions  of  separated  flow  in  high  Mach  num¬ 
ber,  low  Reynolds  number  flows.it  is  necessary  to  include  the  pressure 
variation  across  the  boundary  layer  in  the  formulation  of  the  mathematical 
model.  If  the  effect  of  surface  curvature  on  the  development  and  separa¬ 
tion  of  a  boundary  layer  on  a  curved  surface  is  to  be  calculated, then  the 
normal  momentum  equation  must  be  solved  simultaneously  with  the 
conservation  equations  in  the  streamwise  direction.  In  this  new  theoreti¬ 
cal  analysis  we  have  derived  a  model  to  describe  the  viscous-inviscid 
interaction  for  attached  and  separated  flows  over  two-dimensional  curved 
compression  surfaces  in  hypersonic  flow.  Here  we  have  obtained 
simultaneous  solutions  to  the  integral  forms  of  the  boundary  layer  equations 
for  the  conservation  of  mass,  streamwise  momentum,  moment  of  stream- 
wise  momentum,  normal  momentum  and  energy.  Because  the  formulation 
of  the  analysis  was  an  order  of  magnitude  more  difficult  than  the  earlier 
analysis, we  had  to  use  both  Cohen  and  Reshotko  solutions  and  quartic 
polynomials  to  describe  the  velocity  and  enthalpy  profiles.  We  found 
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that  initial  response  of  the  boundary  layer  to  a  streamwise  pressure 
gradient  at  the  beginning  of  the  interaction  region.  Small  values  of  the 
parameter  ^  can  change  the  response  from  supercritical  to  sub- 

r W 

critical.  Comparisons  among  the  new  analysis,  the  earlier  analysis  where 
d-p/ dy  -  0,  and  the  experimental  measurements  showed  that  introducing  the 
normal  pressure  gradient  gave  larger  streamwise  pressure  gradients  in 
the  separation  and  reattachment  regions,  which  was  in  closer  agreement 
with  experiment.  We  did  not  encounter  any  problems  with  the  integration 
in  the  neck  region,  which  was  the  area  where  the  theory  with  dp/dy  *  0 
broke  down  in  the  strong  interaction  regime.  The  calculations  for  this 
condition  were  in  good  agreement  with  measurements  from  the  experimental 
program.  The  calculations  indicated  that  pressure  ratios  of 

order  of  1.  6  are  generated  in  the  separation  and  reattachment  regions  of 
the  flow. 

The  Experimental  Study  of  Shock  Wave-Boundary  Layer  Interaction  on 
Curved  Compression  Surfaces 

In  this  study,  skin  friction,  heat  transfer  and  pressure  measure¬ 
ments  were  made  on  a  series  of  flat  plate-cylindrical  arc-wedge  compres¬ 
sion  surfaces  to  examine  the  effects  of  radius  of  curvature  on  the 
occurrence  and  structure  of  separated  regions.  Pitot  and  cone  pressure 
traverses  were  also  made  to  obtain  an  estimate  of  the  pressure  gradients 
normal  to  the  surface  in  the  separation  and  reattachment  regions.  In 
these  flows,  the  first  indication  that  flow  separation  was  to  occur  was 
provided  by  the  skin  friction  measurements  on  the  compression  surface; 
the  skin  friction  was  first  observed  to  be  negative  on  the  surface  down¬ 
stream  of  the  flat  plate.  On  the  compression  surfaces  with  large  radii  of 
curvature,  the  incipient  separation  condition  could  not  be  detected  with  any 
accuracy  by  observing  the  first  occurrence  of  an  inflexion  point  in  the 
pressure  distribution  or  a  change  in  the  form  of  the  heat  transfer  distri¬ 
bution,  The  static  pressure  distributions,  evaluated  from  the  pitot  and 
cone  pressure  measurements  across  the  boundary  layers  in  the  separation 
and  reattachment  regions,  indicated  that  the  normal  pressure  gradients 
were  of  the  same  magnitude  as  those  in  the  streamwise  direction.  Although 


IX 


the  experiments  were  conducted  for  a  range  of  ratios  of  60/R  ,  which  were 
similar  to  those  used  by  Kuehn,  «.nd  Sterrett  and  Emery,  we  found  that 
incipient  separation  was  far  less  sensitive  to  surface  curvature  than 
indicated  in  these  earlier  studies.  This  we  believe  results  from  the 
structure  of  the  laminar  boundary  layer  in  hypersonic  flow  over  highly 
cooled  walls  where  the  streamwise  distances  for  a  surface  disturbance 
to  propagate  to  the  edge  of  th<>  boundary  layer  (where  most  of  the  momen¬ 
tum  is  concentrated)  are  greater  than  the  length  of  the  curved  surface. 

_  yfc/ 

The  measurements  were  correlated  in  the  form  0  .  versus  t  (1+ 
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1 .  INTRODUCTION 


The  study  of  the  non- similar  development  of  a  boundary  layer  in  a 
strong  favorable  or  adverse  pressure  gradient  has  provided  one  of  the  more 
challenging  areas  of  research  since  the  formulation  of  the  boundary  layer 
theory.  While  most  theoretical  methods  for  predicting  the  structure  and 
development  of  boundary  layers  are  based  on  the  assumption  that  the  veloc¬ 
ity  profiles  develop  in  a  similar  or  locally  similar  fashion,  the  development 
of  boundary  layers  in  separated  flows  over  bodies  or  in  their  near  wake  can 
be  described  only  by  methods  which  model  the  non-similar  nature  of  the  flow 
mechanics  inherent  in  such  regions.  Within  the  past  decade  this  subject  has 
received  considerable  attention  because  of  its  importance  to  the  design  of 
compression  surfaces  on  high  speed  vehicles  and  the  characteristics  of  the 
wake  of  re-entry  vehicles.  Although  many  experimental  and  theoretical  studies 
of  laminar  and  turbulent  separated  regions  have  been  conducted,  the  under¬ 
standing  of  the  basic  fluid  mechanics  is  far  from  complete.  However,  in  one 
area,  that  of  laminar  separated  flows  caused  by  shock  wave-boundary  layer 
interaction  in  supersonic  flow,  theoretical  analyses  based  on  solutions  to  the 
integral  forms  of  the  conventional  boundary  layer  equations  have  been  success¬ 
ful  in  predicting  many  of  the  characteristics  observed  in  experimental  studies 
of  these  regions. 

In  most  practical  situations  regions  of  shock  wave -boundary  layer 
interaction  are  turbulent  in  supersonic  flight;  however,  in  hypersonic  flows  at 
high  altitudes,  fully  laminar  separated  regions  exist  in  continuum  flow.  In  the 
high  Mach  number,  low  Reynolds  number  ’gime,  where  laminar  regions  of 
shock  wave-boundary  layer  interaction  are  of  practical  importance,  some  of 
the  basic  assumptions  made  in  the  theoretical  methods  successful  in  predicting 
the  characteristics  of  such  regions  in  supersonic  flow  are  invalid.  In  this 
section  we  will  discuss  some  fundamental  features  of  these  theoretical  models 
and  where  they  can  be  changed  to  describe  regions  of  shock  wave-laminar 
boundary  layer  interaction  in  high  Mach  number,  low  Reynolds  number  flows 
over  highly  cooled  curved  compression  surfaces. 
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Within  the  past  few  years  theoretical  methods  based  on  solutions  to 

the  integral  forms  of  the  boundary  layer  equations*"*®  and  finite  difference 
11-15 

solutions  of  the  boundary  layer  equations  have  been  successful  in  pre¬ 

dicting  the  occurrence  and  properties  of  regions  of  laminar  attached  and 
separated  flows  over  adiabatic  and  cooled  compression  and  expansion  corners. 
In  these  analyses  the  common  fundamental  assumption  is  that  the  streamwise 
pressure  gradient  in  such  regions  of  shock  wave-boundary  layer  interaction 
is  generated  solely  by  a  strong  interaction  between  the  viscous  and  inviscid 
flow  fields.  Thus,  the  initial  development  of  this  viscous  layer  must  be  inde¬ 
pendent  of  the  disturbance  which  promoted  the  interaction  whether  it  be  an 
incident  shock  or  compression  or  expansion  corner.  This  "free  interaction" 
assumption  has  been  verified  experimentally  for  small  separated  regions  in 
supersonic  anc  hypersonic  flows.  In  fact,  experiment  has 

shown  that  for  small  separated  regions,  the  distribution  of  heat  transfer  and 
pressure  on  the  wall  over  the  complete  interaction  region  for  wedge  and  exter¬ 
nally  generated  shock  induced  interactions  of  equal  strengths  are  almost 
identical.  This  result  is  of  importance  since  it  verifies  to  some  extent  the 
assumption  made  in  the  theoretical  analyses  that  the  interaction  can  be  treated 
in  terms  of  a  two-point  boundary  value  problem.  Thus,  the  disturbance  which 
promoted  the  interaction  enters  the  problem  via  the  downstream  condition,  and 
a  unique  solution  is  obtained  by  an  iterative  procedure  to  satisfy  both  the  up¬ 
stream  and  downstream  boundary  conditions  simultaneously. 

Most  of  the  analyses  of  shock  wave -boundary  layer  interaction  have 
been  based  on  the  classical  assumption  that  the  static  pressure  remains  con¬ 
stant  across  the  viscous  layer.  Although  this  assumption  can  be  justified  for 
thin  boundary  layers  in  supersonic  flow,  in  the  high  Mach  number,  low  Reynolds 
number  regime  where  laminar  flows  are  of  practical  importance  large  pressure 
differences  can  be  developed  across  a  thick  boundary  layer  in  the  separation 
and  reattachment  regions.  In  these  flows,  as  shown  in  the  present  experimen¬ 
tal  study,  the  length  of  the  separation  region  is  of  the  same  magnitude  as  the 
boundary  layer  thickness,  and  since  large  wedge  angles  are  required  to 
promote  flow  separation,  there  is  pronounced  curvature  of  the  viscous  layer 
in  the  reattachment  zone  (see  Fig.  1).  Under  these  conditions  there  is  a 
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serious  question  as  to  the  validity  of  the  assumptions  made  in  obtaining  the 
conventional  boundary  layer  equations  from  the  Navier-Stokes  equations. 

A  pressure  gradient  is  developed  normal  to  the  surface  as  a  result  of 
curvature  of  the  streamlines  in  the  boundary  layer.  This  curvature  can  be 
generated  by  the  response  of  the  boundary  layer  to  a  longitudinal  pressure 
gradient  induced  by  a  disturbance  downstream,  or  from  the  curvature  of  the 
surface  upon  which  the  boundary  layer  is  growing.  In  general,  the  curvature 
at  the  outer  edge  of  the  boundary  layer  is  not  the  same  as  the  wall  curvature; 
thus,  the  normal  pressure  gradient  is  caused  by  a  combination  of  these  two 
effects.  I’  these  flows  it  is  possible  to  calculate  the  effect  of  surface  curva¬ 
ture  on  boundary  layer  development  and  separation  only  if  a  simultaneous 
solution  is  obtained  to  the  equations  including  the  conservation  equation  in  a 
direction  normal  to  the  surface. 

Similar  solutions  to  the  boundary  layer  equations  including  the  normal 
momentum  equation  have  been  derived  for  a  certain  class  of  surfaces  for 

which  the  boundary  layer  equations  reduce  to  a  set  of  ordinary  differential  equa¬ 
tions.  The  solutions  indicate  that  increasing  surface  curvature  delays  boundary 
layer  separation.  However,  in  these  methods  there  is  no  mechanism  provided 
for  the  interaction  between  the  viscous  and  inviscid  flow  fields,  and  it  is  this 
mechanism  which  has  been  used  in  the  integral  methods  to  describe  the  non- 
similar  development  leading  to  boundary  layer  separation.  In  the  similar 
solution  the  outer  edge  of  the  boundary  layer  assumes  the  curvature  of  the 
surface  and  this  is,  of  course,  unrealistic  for  small  radii  of  curvature. 

One  phenomenon  which  has  arisen  when  the  conventional  boundary 
layer  equations  have  been  used  to  describe  the  viscous-inviscid  interaction  in 
high  h  ach  number  flows  over  highly  cooled  walls  is  the  supercritic  ail  boundary 
layer.  The  concept  of  subcritical  and  supercritical  boundary  layers  was  first 
introduced  by  Crocco  to  describe  the  response-of  the  boundary  layer  when 
subjected  to  an  adverse  pressure  gradient  acting  uniformly  across  the  viscous 
layer.  A  subcritical  flow  was  defined  as  one  in  which  aLS/cLp  was  positive, 
while  this  quantity  was  negative  for  a  supercritical  boundary  layer.  Thus,  a 
subcritical  boundary  layer  will  thicken  when  subjected  to  an  adverse  pressure 
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gradient,  while  a  supercritical  layer  will  thin  under  these  conditions.  Physi¬ 
cally.  a  subcritical  boundary  layer  is  one  in  which  the  increase  in  the  area  of 
the  streamtube  beneath  the  sonic  line  is  greater  than  the  area  contraction  in 
the  supersonic  portion  of  the  layer.  Therefore  a  subcritical  boundary  layer 
can.  by  interacting  with  the  inviscid  stream,  develop  its  own  positive  pressure 
gradient,  whereas  no  such  mechanism  exists  if  the  boundary  layer  is  super- 
c  itical.  Thus,  within  the  conventional  boundary  layer  formulation  there  is 
no  mechanism  by  which  a  supercritical  boundary  layer  can  proceed  smoothly 
to  separation,  a  point  at  which  the  viscous  layer  must  be  subcritical  regard¬ 
less  of  the  local  Mach  number  or  degree  of  wall  cooling. 

In  order  to  overcome  this  problem  within  the  framework  of  conven¬ 
tional  boundary  layer  theory,  it  has  been  suggested  that  the  boundary  layer 
undergo  a  rapid  change  or  jump  to  bring  it  from  a  supercritical  to  a  subcritical 
response  ahead  of  separation.  Grange,  Klineberg,  and  Lees  have  developed 

expressions  equivalent  to  the  normal  shock  relationship  for  the  magnitude  of 
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the  jump,  while  Holden  has  assumed  that  the  jump  is  just  sufficient  to  bring 
the  boundary  layer  to  a  point  where  it  exhibits  a  subcritical  response.  Meas¬ 
urements  in  the  separation  region  of  boundary  layers  which  theoretically  should 
exhibit  a  supercritical  response  do  not  indicate  that  the  mechanism  of  boundary 
layer  separation  differs  significantly  from  that  of  a  subcritical  boundary  layer. 
Certainly  there  does  not  exist  a  discrete  phenomenon  in  the  separation  process 
to  which  a  "shock-like"  characteristic  could  be  ascribed.  We  believe  that  the 
supercritical -subcritical  jump  should  be  viewed  as  a  device  which  must  be 
used  because  of  the  limitations  of  the  integral  solution  to  the  conventional 
boundary  layer  equations,  and  does  not  represent  a  part  of  the  mechanism  by 
which  separation  takes  place. 

29 

It  is  clear  from  the  experimental  studies  of  laminar  regions  of  shock 
wave-boundary  layer  over  highly  cooled  compression  surfaces  that  the  boundary 
layers  which  in  conventional  theoretical  framework  are  classed  as  "super¬ 
critical"  do  not  exhibit  a  supercritical  response.  There  is  no  evidence  of  local 
boundary  layer  thinning  or  an  increase  in  the  heat  transfer  and  skin  friction  at 
the  beginning  of  the  interaction  region.  The  measurements  indicate  that  the 
boundary  layer  thickens  upstream  of  separation,  which  implies  a  larger  area 
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change  of  the  streamtubes  in  the  subsonic  region  relative  to  the  decrease  in 
area  of  the  streamtubes  in  the  supersonic  region  of  the  boundary  layer.  This 
in  turn  implies  that  a  pressure  gradient  exists  across  the  boundary  layer.  It 
should  also  be  noted  that  in  these  flows  the  extent  of  the  upstream  influence 
ahead  of  the  separation  point  (from  skin  friction  measurements)  is  of  the  order 
of  the  boundary  layer  thickness.  Thus>  in  these  flows •  the  assumptions  made  in 
conventional  boundary  layer  that  the  normal  pressure  gradient  can  be  neglected 
are  invalid. 

In  order  to  describe  properly  the  viscous-inviscid  interaction  leading 
to  boundary  layer  separation  in  high  Mach  number,  low  Reynolds  number  flow, 
it  is  necessary  to  include  the  pressure  variation  across  the  boundary  layer 
in  the  formulation  of  the  mathematical  model.  If  the  effect  of  surface  curva¬ 
ture  on  the  development  and  separation  of  a  boundary  layer  on  a  curved  surface 
is  to  be  calculated,  then  the  normal  momentum  equation  must  be  solved  simul¬ 
taneously  with  the  conservation  equations  in  the  streamwise  direction.  In  the 
theoretical  analysis,  which  is  described  in  Section  2  of  this  paper,  we  have 
derived  an  analysis  to  describe  the  viscous-inviscid  interaction  for  attached 
and  separated  flows  over  two  dimensional  curved  compression  surfaces  in 
hypersonic  flow.  Here  we  have  obtained  simultaneous  solutions  to  the  integral 
forms  of  the  boundary  layer  equations  for  the  conservation  of  mass,  streamwise 
momentum,  moment  of  streamwise  momentum,  normal  momentum  and  energy. 
The  velocity  and  enthalpy  profiles  in  the  viscous  layer  have  been  represented 
by  the  Cohen  and  Reshotko  analogs  of  the  Falkner-Skan  equations  and  quartic 
polynomials.  It  was  found  that  the  sub- supercritical  boundary  is  strongly 
dependent  on  the  assumptions  made  about  dp/dy  .  Even  small  values  of 
4^9  were  sufficient  to  change  the  response  of  a  boundary  layer  from 
supercritical  to  subcritical.  Comparisons  are  made  between  theoretical 
solutions  and  experimental  measurements  of  skin  friction,  pressure  and  heat 
transfer  from  the  present  and  previous  experimental  studies. 

In  the  experimental  studies  described  in  Section  3  of  this  paper  we 
examined  the  effect  of  model  curvature  and  free  stream  conditions  on  the  occur¬ 
rence  and  properties  of  separated  regions  over  two-dimensional  curved  surfaces. 
In  this  work,  skin  friction,  heat  transfer  and  pressure  measurements  were 
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made  in  attached  and  separated  regions  on  a  series  of  flat  plate-curved  surface 
vedge  configurations.  The  length  and  radius  of  curvature  of  the  curved  sur¬ 
face  were  varied,  as  well  as  the  length  of  the  flat  plate  and  the  unit  Reynolds 
number  of  the  free  stream.  The  skin  friction  measurements  were  used  to 
detect  the  conditions  and  position  at  which  incipient  separation  occurs;  the 
pressure  and  heat  transfer  distributions  were  examined  at  these  conditions. 
Pitot  and  cone  pressure  measurements  were  made  on  the  separation  and 
reattachment  regions  so  that  an  estimate  of  the  static  pressure  difference 
across  the  boundary  layer  could  be  made.  The  incipient  separation  conditions 
were  correlated  with  model  and  free  stream  parameters. 


2.  THEORETICAL  ANALYSIS 

2.  1  THE  BASIC  EQUATIONS 

In  this  formulation  the  integral  forms  of  the  boundary  layer  equations 
are  used  to  describe  the  non- similar  development  of  attached  and  separated 
laminar  boundary  layers  on  highly  cooled  curved  compression  surfaces  in 
hypersonic  flow.  A  simultaneous  solution  is  obtained  to  integral  forms  of  the 
equations  for  the  conservation  of  mass,  streamwise  momentum,  normal 
momentum,  moment  of  streamwise  momentum  and  energy.  Unlike  the  usual 
boundary  layer  solutions  we  do  not  assume  that  the  pressure  is  constant  across 
the  boundary  layer. 

The  boundary  layer  equations  for  two-dimensional  flow  with  a  surface 
oriented  coordinate  system  are  used,  as  shown  in  Fig.  2,  retaining  the 
important  higher  order  terms 

JL  (pU)t  -  (pr)  =  0 

dz  dy 

continuity  (1) 

----If' ^#9 

%-  momentum  (2) 
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energy  (4) 


Here  we  have  assumed  that  since  most  of  the  momentum  is  contained 
in  the  outer  region  of  the  viscous  layer  in  a  hypersonic  highly  cooled  flow,  we 
can  neglect  the  contribution  of  viscous  terms  to  the  normal  pressure  gradient. 
In  addition,  the  assumption  has  been  made  that  the  radius  of  curvature  of  the 
surface  is  large  compared  with  the  thickness  of  the  boundary  layer. 


Integrating  Equation  1  across  the  boundary  layer  and  defining 
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and  assuming  the  isentropic  relationships  for  the  flow  at  the  outer  edge  of  the 
viscous  layer,  Equation  5  becomes 


(1  *  ”?«)  tan  Q 

»re  (i  '  "too)  S* 


d 

dX 


t  iL  *  2 


0*swt) 


t 


d // 


*£jl 

dx 


dx 


('*.*)[*» -(jrj) 


(1  +  *Te) 


Jr-i 

r-i 


+ 


»’e 

1+  me 


H 


— ~  ’  w 

"led  +  ife)  _ 


q  d  t 
"V  dX 


(6) 


7 


where 


Transforming  Equation  2  and  integrating  across  the  boundary  layer, 
we  obtain 


Transforming  Equation  2,  multiplying  through  by  U  ,  then  integrating  across 
the  boundary  layer,  we  obtain 
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and  the  definitions  of  the  velocity,  enthalpy  and  pressure  profile  parameters 
a.  ,  b  and  c  are  given  in  the  Section  2.2.  Integrating  the  normal  moment 
equation  across  the  boundary  layer, 


where 


■jr(—  -  — )<*  Y  =  *r"\1  •  W~  7  S*  -  w /  A 1 

^eW,  Uj  i  \u-Trr,)  ]  ct  X.  d  X, 

0 

~  b's>' -  £  7T  7V~\P  +  (*rh'M')p»rpr> 

-p<"  p.  {£<**-*  7T £r]p] 


tcW -f 


where 
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P"'  $  I.  f‘ *1  dr!  '  •“> /I  ‘ 

P>5  '  ^  J  *1  '>  -  -£T  f  0~  f  <*rl  i  PQ  -f(,-  f')4rj  • 

The  energy  equation  (4),  when  integrated  across  the  boundary  layer,  becomes 


T*.  s.drtT*.  .  St  (>1.)  0) 

TJX1°>'S‘  +  T  tx  /oS'  M'  a.M'S;  S„  OV"U/3/“.  / 


where 


U  Y 


6:  j  -w  ^ 

Equations  (b),  (7),  (8),  and  (10)  c  n  be  rewritten  in  the  form 
<=.,  £  s:  - c 


%  Me  rC 


eta 


^  j.  r  ^ C  -  T)  ■ 

3/ax  '  5f  “  /  ' 


y-  c „ 


Thus,  we  obtain  a  system  of  five  simultaneous  non-linear  differential  equations 

where  c  and  -CL  are  functions  of  the  profile  parameters  a.  ,  b  ,  c  and  M, 

/  /  r  s 

and  5*  .  The  term  /  Jt.  iXl  .  J/_)  d  v  contains  second  derivatives  with 

Jo  ue  '  Ut  Ut' 

respect  to  X  .  However,  since  its  magnitude  is  small,  it  was  calculated  by 
taking  the  forward  difference  of  the  first  derivative. 


L.l  VELOCITY  AND  ENTHALPY  PROFILES 

In  the  earliest  attempts  to  apply  momentum  integral  techniques  to 

describe  the  non- similar  development  of  boundary  layer  in  strong  pressure 

gradients,  the  distributions  of  velocity  and  temperature  were  represented  by 

polynomials  in  terms  of  the  non-dimensionalized  distance  from  the  wall  and 

30 

a  shape  factor.  Subsequently,  Bray,  Gadd  and  Woodger,  and  Lees  and 

Reeves^  showed  that  for  separated  regions  induced  by  shock  wave -boundary 

layer  interaction  on  an  adiabatic  wall,  the  compressible  flow  analogs  of  the 

3  1 

Falkner-Skan  profiles  obtained  by  Cohen  and  Reshotko  better  represented 
ihe  form  of  the  velocity  profiles.  This  approach  was  extended  successfully  to 
flows  with  heat  transfer  by  Holden,^  and  later  by  Klineberg  and  Lees,  ^  by  using 
both  the  velocity  and  enthalpy  profiles  from  the  Cohen  and  Reshotko  solutions. 
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In  the  present  analysis  the  Cohen  and  Reshotko  profiles  have  been 
used  to  describe  the  distribution  of  velocity  and  enthalpy  in  the  attached  and 
separated  viscous  layer.  In  addition,  because  of  the  complexity  of  the  present 
problem,  we  introduce  quartic  polynomials  to  represent  the  form  of  the  veloc¬ 
ity  and  total  enthalpy  profiles  in  the  terms  which  arise  when  the  pressure  is 
allowed  to  vary  across  the  boundary  layer.  The  equations  have  been  formu¬ 
lated  in  such  a  way  that  the  terms  in  which  a  pressure  distribution  must  be 
specified  represent  the  difference  between  the  analysis  for  dp/dy-0  and 

t  0  .  The  quartic  distributions  for  the  velocity  and  enthalpy  profiles 
used  in  these  terms  are  matched  to  the  Cohen  and  Reshotko  profiles  used  in 
remaining  and  more  important  terms.  The  polynomial  forms  chosen  for  the 
velocity  and  enthalpy  profiles  were 

=  1*  (*>  ~  &  1  *  ^ r\2)  *  A  7  (/- rj)’ 

and 

—  =  1  -  'J  ((>-  +  3  *1  )  ~  B  rj  (1-  fj ) 

Here  A  and  8  are  the  slopes  of  the  velocity  and  total  enthalpy  profiles  at  the 

wall,  respectively.  This  particular  representation  was  chosen  so  that  A (X  ) 

and  8  (Xl  could  be  related  directly  to  the  profile  parameters  a,(X)  and  b(x) 

which  are  used  to  characterize  the  velocity  and  enthalpy  profiles  obtained  from 

the  Cohen  and  Reshotko  solutions.  Here  a,(X)  ,  the  velocity  profile  parameter, 

is  defined  as  &&/(**)  for  a^ached  flow  and  Y YS.  for  separated  flow,  where 
„  *(Y/6i ) 

Y  is  the  distance  between  the  wall  and  the  point  of  zone  velocity  in  the 
separated  shear  layer.  b(X)is  defined  as  the  slope  of  the  total  enthalpy  pro¬ 
file  at  the  wall  in  both  attached  and  separated  flows.  The  Cohen  and  Reshotko 
and  polynomial  velocity  profiles  are  linked  by  the  relationship  A  -  a  for 
attached  flow,  and  A  =  &((>-  8a,  *  3  a  1 )/ (1  -  <*•/  for  separated  flow,  (match¬ 

ing  the  point  of  zero  velocity  in  the  separated  shear  layer).  For  the  enthalpy 
profiles  8  =  8  rj  * ^  in  both  attached  and  separated  flows,  where  >]**  is  the 

value  of  rj  at  'Jj Ue  =0.99-  A  comparison  between  the  Cohen  and  Reshotko 
and  polynomial  distributions  is  shown  in  Fig.  3.  Here  it  can  be  seen  that  the 
distributions  are  in  close  agreement  for  attached  boundary  layers,  but  show 
noticeable  departure  with  significant  reverse  flow.  It  is  of  interest  to  compare 
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the  form  of  the  polynomial  and  Cohen  and  Reshotko  velocity  distributions  in 
separated  regions  for  identical  velocity  gradients  at  the  wall,  also  shown  in 
Fig.  3.  The  significant  difference  between  the  profiles  correlated  in  this  way 
is  responsible  in  part  for  discrepancies  between  the  solutions  obtained  from 
integral  methods  using  single  parameter  polynomials  and  Cohen  and  Reshotko 
profiles . 

The  pressure  distribution  across  the  boundary  layer  is  assumed  to  be 
of  the  form 

-  1  -  c  (l  -  r}1  ) 

P  ' 

Pw  -  p . 

where  c  =  — — -  .  This  form  satisfies  the  boundary  conditions  at  the  wall 

and  edge  of  the  boundary  layer,  the  momentum  equation  at  the  wall,  and  is  in 

agreement  with  the  form  of  the  static  pressure  measurements  of  Stroud  and 
3  2 

Miller  across  boundary  layers  on  curved  surfaces.  The  density  distribution 

is  calculated  from  the  equation  of  state  -£■  =  where  -4--  (1 1  S){l  *■  »7e)-7^ef' 

'  e  re  **■ 

Because  of  the  algebraic  complexity  of  evaluating  many  of  the  integrals 
associated  with  the  normal  pressure  gradient,  both  the  manipulation  and  exact 
integration  were  performed  on  a  digital  computer  using  Formax.  These  inte¬ 
grals  were  expressed  exactly  in  terms  of  A  ,  B  ,  C  ,  ->rje  and  5  w  . 

The  integral  quantities  H  ,  J  ,  P  ,  R  ,  etc.,  which  were  determined 
from  the  Cohen  and  Reshotko  solutions  were  fitted  as  a  power  series  in  a-  by 
the  method  of  least  squares  and  expressed  in  the  form 

H  -  cD  i-  c,a.  +■  clai  .  .  .  . 

Similarly,  the  integrals  which  depended  on  the  enthalpy  profiles  were  fitted  as 
a  power  series  in  b  .  Functions  of  both  a.  and  b  were  written  in  terms  of 
component  expressions  composed  of  power  series  in  a.  and  b  . 
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2.3 


SOLUTION  OF  THE  DIFFERENTIAL  EQUATIONS  AND  COMPARISON 
WITH  EXPERIMENT 


The  equations  for  the  continuity  of  mass,  x-  momentum,  moment  of 
X-  momentum,  y- momentum  and  energy  (Equations  6,  7,  8,  9  and  10)  can  be 
written  as  a  system  of  five  simultaneous,  non-linear,  ordinary  differential 
equations  of  the  form 


dh 

dx 


i-  C. 


dc 

dX 


=  Z> , 


where  the  displacement  thickness  o L  ,  the  local  Mach  number  /Vc  and  the 
velocity,  enthalpy  and  pressure  profile  parameters  ^  ,  b  and  c  are  the 
independent  variables  for  which  a  solution  must  be  obtained  at  each  X  station. 
The  algebraic  expressions  for  c  ,  c2  .  ,  etc.,  are  contained  in  Appendix  1 


These  equations  were  solved  within  the  framework  of  a  two-point 

boundary  value  problem  (as  discussed  in  the  Introduction)  where  certain  of  the 

boundary  conditions  upstream  and  downstream  of  the  corner  interaction  are 

known.  A  unique  solution  is  obtained  by  iterating  a  trial  integration  until  all 

the  boundary  conditions,  both  upstream  and  downstream  of  the  interaction 

region,  are  satisfied  simultaneously.  The  integration  was  performed  on  an 

IBM  360  computer  using  a  modification  of  the  5th  order  Runge-Kutta-Merson 

33 

method  developed  at  CAL  by  Partee.  A  complete  integration  from  the 
upstream  to  the  downstream  conditions  took  less  than  one  minute,  with  a 
typical  unique  solution  requiring  approximately  15  iterations. 


The  integration  is  initialized  just  upstream  of  the  corner  interaction 
region.  Here  we  assume  the  development  of  the  velocity  and  enthalpy  profiles 
to  be  self-similar,  hence  -  0  .  If  the  beginning  of  the  interaction 

occurs  over  the  flat  plate  and  the  interaction  between  the  viscous  and  inviscid 
flow  in  this  region  is  "weak"  (  L  small),  then  yyy  =  0  .  Under  this  condi¬ 

tion,  a  solution  for  the  initial  conditions  yields  the  Blasius  solution  for  the 
velocity  and  enthalpy  distribution  in  the  boundary  layer.  If  the  interaction 
over  the  flat  plate  is  "strong"  or  the  boundary  layer  is  subjected  to  a  pressure 
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gradient  other  than  that  associated  with  the  corner  interaction,  then  for 

A  oc  z"  , 

Pm 


(  2  ) 

I  X  /-/ 

\  7*7  *■  l) 

Me  dX  '  X  1 

U  /-  /  / 

where  U e  -  c  X  m  .  Hence  the  velocity  and  enthalpy  profiles  at  the  begin¬ 
ning  of  the  interaction  can  be  related  to  one  of  the  family  of  self-similar 
profiles  obtained  from  the  Cohen  and  Reshotko  solutions. 

In  order  to  proceed  with  a  continuous  integration  toward  the  separation 

c£  cl  cti 

point,  the  initial  response  of  the  boundary  must  be  such  that  ^  and 

he  negative  while  and  are  positive,  i.e.,  the  boundary 

layer  responds  in  a  subcritical  manner.  The  subject  of  subcritical  and  super¬ 
critical  boundary  layers  has  been  discussed  at  some  length  in  the  Introduction. 
We  noted  the  classification  was  originally  derived  from  analyses  in  which  the 
classical  assumption  that  d-p/3y-0  was  made,  and  suggested  that  in  high 
Mach  number,  low  Reyi  olds  number  flows  over  highly  cooled  walls  a  pres¬ 
sure  gradient  generated  across  the  boundary  layer  plays  an  important  role  in 
the  mechanism  of  flow  separation.  In  the  present  analysis  we  can  not  only 
•  nlculatc  the  pressure  difference  developed  across  the  boundary  layer,  but 
also  investigate  the  effect  of  a  non-zero  pressure  gradient  across  the  boundary 
layer  on  its  initial  response.  We  have  made  calculations  of  the  effect  of  small 
positive  values  of  c  on  the  superc  ritical -  subc  ritical  boundary  and  found  that 
it  is  extremely  sensitive  to  tins  parameter.  For  one  of  the  examples  studied 
in  earlier  work  (Rel.  29),  the  separation  of  a  laminar  boundary  layer  in  Mach  20 
airflow  in  the  strong  interaction  regime,  a  value  of  C  of  just  less  than  0.04 
changed  the  response  of  the  boundary  from  supercritical  to  subcritical.  We 
propose  that,  rather  than  using  a  jump  condition  to  overcome  the  problems 
associated  with  the  occurrence  of  a  supercritical  boundary  layer,  a  more  sa*is- 
lactory  approach  from  both  the  physical  and  mathematical  viewpoints  is  to 
specify  a  pressure  difference  across  the  boundary  layer  of  just  sufficient 
magnitude  to  evoke  a  subcritical  response.  Since  the  initial  value  for  C  is 
iar  smaller  than  those  subsequently  calculated  through  the  interaction  region, 
its  only  effect  is  to  perturb  the  boundary  layer,  causing  the  subsequent 
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subcritical  development.  At  the  downstream  end  of  the  interaction,  we 
require  both  a  and  b  to  approach  the  Blasius  profiles  as  the  Mach  number 
approaches  the  prescribed  value  on  the  wedge. 

A  typical  solution  for  the  basic  variables  a.  ,  b  ,  c  and  Me  is  shown 
in  Fig.  4.  Here  the  calculations  were  performed  for  comparison  with  pres¬ 
sure  and  heat  transfer  measurements  made  in  a  Mach  10  airflow  over  a  highly 
cooled  flat  plate-wedge  configuration,  shown  in  Fig.  5.  The  pressure,  heat 
transfer  and  skin  friction  are  related  to  the  basic  variables  by  the  relation¬ 
ships 


At  the  beginning  of  the  interaction  the  Mach  number  is  less  than  the  free  stream 
value  because  of  the  weak  interaction  at  the  leading  edge,  and  consequently 
the  pressure  pg  at  this  point  is  slightly  larger  than  the  static  pressure  of  the 
free  stream.  The  Mach  number  falls  rapidly  with  the  gradient  reaching  a 
maximum  at  separation,  which  is  reflected  in  an  inflexion  in  the  pressure 
distribution  at  this  point.  Both  a  and  b  ,  and  hence  and  ,  decrease 

monatonically  as  separation  occurs,  with  a  and,  of  course,  Cf  ,  equal  to 
zero  at  separation.  In  the  separated  region,  a.  (now  re-defined  as  Y*/6*  ) 
increases  to  a  maximum  at  the  point  of  maximum  displacement,  the  point 
where  b  and  the  heat  transfer  rate  reach  a  minimum.  In  this  region  the  Mach 
number  approaches  a  constant  value,  causing  a  region  of  approximately 
constant  pressure  --  the  plateau  pressure  region.  Downstream  of  the  plateau 
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region  the  Mach  number  decreases  through  the  reattachment  compression 
process  to  approach  asymptotically  the  wedge  Mach  number,  a-  returns  to 
zero  at  the  reattachment  point,  and  then  in  the  attached  flow  downstream  CL 
and  b  increase  to  approach  their  Blasius  values  on  the  wedge.  We  note  that 
reattachment  occurs  approximately  halfway  up  the  reattachment  compression 
rise  and  the  heat  transfer  is  a  maximum  at  the  end  of  the  compression  proc¬ 
ess.  At  the  beginning  of  the  interaction  we  assume  that  c  is  zero  or  very 
small.  C  increases  through  the  separation  process,  reaching  a  local  maxi¬ 
mum  at  the  point  of  boundary  layer  separation.  From  Fig.  5  it  can  be  seen 
that  the  effect  of  including  the  normal  pressure  gradi  »nt  is  to  accentuate  the 
pressure  gradient  in  the  separation  region.  Because  C  ,  and  hence  the  pres¬ 
sure  difference  across  the  boundary  layer,  drops  rapidly  in  the  separated 
region,  the  new  theory  gives  a  more  well-defined  plateau  pressure  region  than 
the  theory  with  dp/dy-0  .  Downstream  of  the  maximum  displacement  thick¬ 
ness  c.  increases  to  reach  a  local  maximum  at  the  reattachment  point.  Again 
the  reattachment  compression  rise  is  steeper  for  the  theory  including  dp/ dy  . 
c  approaches  zero  as  the  curvature  of  the  shear  layer  decreases.  The  pres¬ 
ent  theory  predicts  slightly  larger  heat  transfer  rates  in  the  separated  region 
than  the  measurements  or  the  theory  with  dp/dy  =  0  indicate;  however,  it  is 
in  better  agreement  with  the  measurements  downstream  of  reattachment  than 
the  more  simple  theory. 

In  the  above  comparison  we  have  seen  that  the  effect  of  including  the 
normal  pressure  gradient  in  the  theory  is  to  steepen  the  gradient  in  the  sepa¬ 
ration  and  reattachment  regions,  giving  a  better  defined  plateau  region. 
However,  since  the  boundary  layer  was  comparatively  thin,  the  maximum 
pressure  ratio  across  the  boundary  was  less  than  1.25.  In  the  experimental 
study  presented  in  this  paper,  the  boundary  layers  were  an  order  of  magnitude 
thicker,  the  Mach  number  was  a  factor  of  1.5  larger,  and  as  a  consequence, 

the  pressure  rise  required  to  separate  the  boundary  layer  and  the  associated 
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turning  angles  of  the  viscous  layer  were  significantly  increased.  Earlier, 
when  we  attempted  to  make  a  comparison  between  similar  experimental  meas¬ 
urements  and  the  theory  without  allowing  for  a  normal  pressure  gradient, 
we  found  we  were  unable  to  integrate  completely  to  the  end  of  the  interaction 
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region  and  thus  obtain  a  unique  solution.  This  problem  was  associated  with 
the  derivatives  becoming  discontinuous  just  downstream  o£  reattachment,  and 
at  this  time  it  was  suggested  that  a  possible  cause  of  this  was  the  failure  to 
describe  the  flow  accurately  by  allowing  for  normal  pressure  gradients.  This 
was  a  stimulus  toward  developing  the  theory  with  dp/dy  £  0  .  Using  the 

present  analysis,  we  were  able  to  obtain  complete  solutions  where  the  earlier 
theory  had  failed.  A  comparison  between  the  theory  and  the  measurements  of 
skin  friction,  pressure  and  heat  transfer  made  on  the  flat  plate  -  18*  wedge  in 
a  Mach  14.  5  airflow  is  shown  in  Fig.  6.  The  calculated  distribution  of  the 
static  pressure  ratio,  shown  in  Fig.  6a.  indicates  that,  as  observed  earlier, 
this  ratio  is  greatest  at  the  separation  and  reattachment  points;  however,  in 
this  case,  pressure  ratios  greater  than  1.6  were  generated  across  the  boundary 
layer.  Since  the  length  of  the  separated  region  is  comparatively  small,  a 
pressure  gradient  still  exists  across  the  viscous  layer  at  the  point  of  maximum 
displacement  thickness.  The  predicted  skin  friction  distribution  is  in  agree¬ 
ment  with  the  measurements,  though  the  theory  underestimates  the  length  of 
the  separated  region.  The  maximum  negative  value  occurs  on  the  wedge  as 
the  experiments  indicate,  and  the  predicted  maximum  occurs  at  the  end  of  the 
reattachment  compression  process  and  is  in  close  agreement  with  the  meas¬ 
urements.  The  comparison  between  the  measured  pressure  distribution  and 
the  predictions  for  the  pressure  at  the  edge  of  the  boundary  layer  and  on  the 
wall,  given  in  Fig.  6b,  show  clearly  the  importance  of  the  variation  in  pres¬ 
sure  across  the  boundary  layer  in  the  analysis.  Only  by  including  the  pressure 
rise  across  the  boundary  layer  can  the  measured  pressure  rise  in  the  separa¬ 
tion  and  reattachment  regions  be  predicted.  As  we  observed  in  the  comparison 
with  the  Mach  10  data,  the  theory  overestimates  the  heat  transfer  rates  in  the 
separated  region  (see  Fig.  6c);  however,  the  level  of  heating  in  the  reattach¬ 
ment  region  is  predicted  to  good  accuracy. 
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3.  EXPERIMENTAL  STUDIES 


3.1  INTRODUCTION 

In  the  experimental  study,  we  exa  mined  the  effect  of  surface  curva¬ 
ture  and  Reynolds  number  on  the  occurrence  and  properties  of  separated 
laminar  boundary  layers  on  highly  cooled  compression  surfaces  in  hypersonic 
airflow.  Here,  in  addition  to  detecting  the  conditions  required  to  promote 
incipient  separation,  by  examining  the  position  of  the  point  of  zero  skin  fric¬ 
tion  and  the  distribution  of  pressure  and  heat  transfer  we  hoped  to  obtain  a  better 

understanding  of  the  mechanism  of  separation  on  curved  compression  surfaces. 

1 6  34 

Earlier  studies  by  Chapman,  Kuehn  and  Larson,  Kuehn,  and  Sterrett  and 
35 

Emery  in  supersonic  flow  had  shown  that  the  pressure  rise  required  to  pro¬ 
mote  incipient  separation  increases  markedly  when  a  curved  surface  is 
introduced  to  replace  the  abrupt  junction  of  a  flat  plate-wedge  compression 
surface.  These  measurements,  which  were  made  in  both  laminar  and  turbu¬ 
lent  boundary  layers  over  adiabatic  surfaces,  indicated  that  the  surface 
curvature  can  increase  the  pressure  rise  to  promote  incipient  separation  by 
as  much  as  a  factor  of  3.5. 

Analytical  studies  of  the  effect  of  surface  curvature  on  boundary  layer 
separation  have  indicated  that  concave  curvature  tends  to  inhibit  flow  separa¬ 
tion  whereas  convex  curvature  produces  the  opposite  effect.  An  important 
parameter  controlling  the  magnitude  of  this  effect  was  found  to  be  the  inverse 
root  of  the  Reynolds  number  based  on  the  local  radius  of  curvature.  These 
analyses  are  restricted  to  a  class  of  similar  surfaces  for  which  the  second- 
order  boundary  layer  equations  reduce  to  a  set  of  ordinary  differential  equations 
for  which  similar  solutions  exist.  In  addition,  it  is  assumed  that  8/R*-<-\ 
and  there  is  no  interaction  between  the  viscous  and  inviscid  flow;  thus,  the  edge 
of  the  boundary  layer  follows  the  surface  curvature  and  there  is  no  mechanism 
for  the  non-similar  development  of  the  boundary  layer  leading  to  separation. 

In  hypersonic  low  Reynolds  number  flow,  the  interaction  between  the 
nor.- similar  growth  of  a  laminar  boundary  layer  and  the  inviscid  flow  around  a 
concave  or  convex  surface  is  of  great  importance.  In  addition  to  the  viscous- 
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inviscid  interaction  which  would  be  induced  in  the  absence  of  surface  curva¬ 
ture,  we  have  additional  transverse  pressure  gradients  generated  as  the 
boundary  layer  attempts  to  follow  the  curvature  of  the  surface.  The  inter¬ 
play  between  the  viscous  and  inertial  forces  generated  by  viscous-inviscid 
interaction  and  surface  curvature  control  the  development  of  the  viscous  layer 
and  the  conditions  under  which  separation  occurs. 

In  the  present  experimental  study,  we  attempted  to  explore  the  condi¬ 
tions  under  which  the  inertial  forces  generated  by  surface  curvature  become 
of  importance.  Beginning  with  a  flat  plate-wedge  compression  surface,  sur¬ 
face  curvature  was  introduced  by  interposing  cylindrical  sections  between  the 
flat  plate  and  wedge.  A  series  of  circular  cylindrical  sections  with  different 
radii  of  curvature  was  tested;  each  cylindrical  section  contained  segments 
which  could  be  added  or  subtracted  to  change  the  degree  of  turning.  The 
position  and  angle  of  the  wedge  was  adjusted  so  that  it  formed  a  continuous 
surface  tangential  to  the  trailing  edge  of  each  cylindrical  segment.  In  the 
following  section,  the  test  facilities,  model  geometry  and  instrumentation  are 
described  and  the  results  of  the  studies  of  the  flow  establishment  times  and 
finite  span  effects  are  reviewed.  The  results  from  the  experimental  measure¬ 
ments  are  then  discussed. 

3.2  EXPERIMENTAL  PROGRAM 

This  experimental  program  was  conducted  in  the  CAL  4-ft  shock 
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tunnel.  This  tunnel,  which  operates  on  tailored  interface  principles,  uses 
a  heated  helium  driver.  For  the  test  conditions  chosen  in  these  studies  the 
useful  running  time  was  approximately  5  milliseconds. 

Photographs  of  the  models  are  shown  in  Fig.  7.  The  models  consisted 
of  a  flat  plate  followed  by  a  circular  cylindrical  surface  which  terminated  into 
a  second  flat  plate.  Cylindrical  surfaces  with  radii  of  curvature  of  6,  12  and 
20  inches  were  inserted  between  the  flat  plate  and  ramp  in  the  program  and 
these  surfaces  were  constructed  in  segments  so  that  the  degree  of  turning 
could  be  varied.  The  second  flat  plate,  or  ramp,  was  supported  so  that  the  slope 
at  the  junction  between  the  curved  surface  and  the  ramp  was  continuous.  The 


19 


length  of  the  leading  edge  of  the  flat  plate  could  be  varied  to  give  lengths  of 
flat  plate  from  12  to  17  inches. 

The  flat  plate,  cylindrical  sections,  and  ramp  were  fully  instrumented 
with  skin  friction,  heat  transfer  and  pressure  transducers  (see  Fig.  7).  The 
thin  film  heat-transfer  gages,  which  were  coated  with  an  insulating  film  of 
titanium  dioxide,  were  used  in  conjunction  with  analog  networks  which  converted 
their  output  directly  into  heat  transfer  rate.  Both  orifice-type  and  flush- 
mounted  pressure  transducers  were  used  in  the  model.  These  exposed 
diaphragm -type  piezoelectric  transducers,  described  in  Ref.  37.  use  lead 
zirconium  titanate  crystals  and  have  a  nominal  sensitivity  of  2  volts/psi.  The 
skin  friction  transducers,  described  in  Ref.  38,  consist  of  a  1/4-inch  diameter 
diaphragm  mounted  flush  with  the  surface  of  the  model,  with  the  diaphragm 
connected  by  a  fixture  to  a  crystal  mounted  in  bending.  This  transducer  which 
is  compensated  for  normal  pressure,  acceleration  and  thermal  effects  has  a 
nominal  sensitivity  of  20  volts/psi. 

The  pitot  and  cone  pressure  probes  used  in  the  rake  were  1/ 8-inch  in 
diameter  and  had  the  diaphragm  of  the  transducer  mounted  directly  behind  the 
orifice  to  ensure  a  maximum  response.  The  transducer  is  compensated  for 
thermal  and  acceleration  loads  and  directly  connected  to  an  F.E.T.  circuit  to 
give  a  low  output  impedance.  A  typical  sensitivity  is  50  mV/psi.  A  double¬ 
pass  schlieren  system  with  a  focal  length  of  10  ft  was  used  to  obtain  photographs 
of  the  flow  over  the  models. 

The  test  program  was  conducted  at  free  stream  Mach  numbers  of  14.0 
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and  14.5  and  Reynolds  numbers  of  7.2  x  10  and  1.8  x  10  per  ft,  respectively. 
Reservoir  pressures  of  1500  and  4200  psia  were  used  at  a  common  stagnation 
temperature  of  5330*R. 

3.3  TIME  ESTABLISHMENT  AND  FINITE  SPAN  EFFECTS 

Two  important  considerations  in  a  shock  tunnel  study  of  two-dimensional 
compression  surfaces  are  (a)  is  the  flow  steady?  and  (b)  are  the  measurements 
made  in  a  "two-dimensional  flow"?  Both  of  these  questions  have  received 
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careful  and  extensive  study  during  the  course  of  the  measurements  in  attached 

and  separated  regions  on  two -dimentional  models.  Since  many  of  the  detailed 
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results  of  these  studies  have  been  published  in  earlier  papers,  *  7  rather 
than  present  similar  measurements  made  in  the  present  study,  we  will  sum* 
marize  the  combined  conclusions. 

For  free  stream  conditions  and  model  configurations  at  which  the 
boundary  layer  remained  completely  attached,  the  viscous  and  inviscid  flows 
were  established  within  the  starting  time  of  the  tunnel  (less  than  1  millisecond). 
The  measurement  indicates  that  it  takes  no  longer  to  establish  an  attached 
boundary  layer  in  a  strong  pressure  gradient  than  it  does  to  establish  the 
boundary  layer  over  a  flat  plate.  The  studies  at  Mach  numbers  from  14  to  20 
and  unit  Reynolds  numbers  from  2, 000/ in  to  14,000/in  indicated  that  separated 
regions  can  take  as  long  as  2.5  milliseconds  to  establish.  Our  measurements 
indicate  that  during  flow  establishment  the  leading  edge  of  the  separated  region 
propagates  forward  at  approximately  the  speed  of  sound  of  the  gas  close  to  the 
Wc\ll.  However,  we  have  found  that  in  some  cases  boundary  layers  just  at  the 
incipient  separation  condition  can  take  almost  as  long  to  stabilize  as  a  large 
separated  region. 

The  effect  of  finite  model  span  on  the  properties  of  both  attached  and 
separated  boundary  layers  was  investigated  by  varying  the  model  span  and 
adding  side  fences  while  monitoring  both  the  longitudinal  and  spanwise  distribu¬ 
tions  on  the  model.  Although  increasing  the  model  span  from  1  ft  to  1.5  ft  did 
not  measurably  influence  the  properties  of  attached  interaction  regions,  it  did 
cause  a  small  increase  in  the  size  of  large  separated  regions.  However,  fur¬ 
ther  increasing  the  span  to  2  ft  had  little  additional  effect.  Adding  side  fences 
to  this  latter  configuration  did  influence  the  distribution  of  properties  on  the 
model.  However,  by  adding  side  fences  to  the  1—  ft  span  model  we  obtained  a 
larger  separated  region  than  that  on  the  2-ft  span  model  with  or  without  side 
fences.  Side  fences  should  be  used  with  care,  for  by  adding  them  it  is  possible 
to  change  the  basic  flow.  We  found  that  a  uniform  distribution  across  the 
center  span  of  the  model  is  not  a  good  criterion  to  determine  whether  the  flow 
is  "two-dimensional."  The  criterion  which  has  been  adopted  here  has  been  to 
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vary  the  transverse  boundary  conditions  by  adding  spanwise  extensions  and 
side  fences  until  successive  change  produced  no  measurable  effect  on  the 
longitudinal  or  transverse  distribution  of  skin  friction*  heat  transfer  or  pres¬ 
sure  on  the  model. 

3.4  SKIN  FRICTION,  HEAT  TRANSFER  AND  PRESSURE  DISTRIBUTIONS 

During  the  experimental  program,  skin  friction,  heat  transfer,  pressure 
and  schlieren  measurements  were  made  on  the  flat  plate-cylindrical  surface- 
ramp  models  with  radii  of  curvature  of  0,  6,  12  and  20  inches  and  ramp  angles 
from  15*  to  35*.  All  the  measurements  made  at  a  particular  configuration 
and  test  condition  were  recorded  simultaneously,  and  in  some  cases  1  15  data 
channels  were  required.  Since  insufficient  space  prohibits  showing  here  all 
of  these  measurements,  series  of  distributions  are  presented  which  are 
representative  of  the  important  characteristics  observed  in  the  study. 

The  distributions  of  properties  typical  of  those  found  in  attached  and 

separated  regions  over  flat  plate-wedge  compression  surfaces  are  shown  in 
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Fig.  8.  Here,  as  in  earlier  studies  ,  the  first  evidence  that  the  boundary 
layer  was  close  to  separation  was  exhibited  in  the  character  of  the  skin  fric¬ 
tion  distribution  downstream  of  the  flat  plate-wedge  junction.  Fig.  8a  shows 
that  the  point  of  minimum  skin  friction  occurs  on  the  ramp  and,  as  found 
earlier  and  observed  in  Fig.  8b,  the  separation  bubble  forms  on  the  ramp 
and  feeds  forward  as  the  wedge  angle  increases.  We  observe  that  in  attached 
flows  the  influence  of  the  wedge  extends  upstream  less  than  one  boundary 
layer  thickness,  and  there  is  no  evidence  to  support  the  speculation  that 
a  small  separated  region  is  always  present  at  the  flat  plate-wedge  junction. 

When  separation  occurs  (as  shown  in  Fig.  8b),  the  initial  development  of 
this  region  is  asymmetric  with  the  small  separation  bubble  displaced  over 
the  ramp.  The  pressure  and  heat  transfer  distributions  are  typical  of  those 
close  to  incipient  separation,  the  pressure  measurements  exhibiting  a  distinct 
inflexion  point,  while  the  heat  transfer  in  the  separation  region  is  reduced 
significantly  below  the  flat  plate  value.  Increasing  the  wedge  angle  to  24® 

(Fig.  8c)  promotes  a  large  region  of  separated  flow  situated  symmetrically 
about  the  hinge  line.  The  pressure  distribution  exhibits  the  classical  constant 
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pressure  region,  and  it  is  interesting  to  note  that  in  this  region  the  skin  friction 
also  remains  constant,  indicating  that  the  structure  of  the  recirculation  region 
consists  of  one  rather  than  a  number  of  counter  rotating  vortices.  In  this 
separated  region  the  pressure  rise  to  reattachment  is  0.58  of  the  maximum 
pressure  rise  from  the  plateau  pressure.  The  heat  transfer  reaches  a  maxi¬ 
mum  at  the  end  of  the  reattachment  compression  process,  and  downstream  of 
thid  point  the  surface  pressure  relaxes  to  the  wedge  pressure  through  an 
expansion  process  generated  when  the  compression  waves  on  the  separation 
and  reattachment  regions  coalesce  to  form  the  wedge  shock. 

Introducing  circular  cylindrical  sections  of  6-inch  radius  between  the 
flat  plate  and  wedge  had  only  a  small  effect  on  the  wedge  angle  to  promote 
incipient  separation  and  the  structure  of  the  separated  regions  on  configura¬ 
tions  which  promoted  separated  flow.  The  distributions  of  surface  properties 
for  two  of  the  configurations  tested  are  shown  in  Figs.  9a  and  9b.  Comparing 
9a  with  8b,  we  see  that  for  the  18*  ramp  angle,  introducing  the  curved  section 
caused  a  slight  reduction  in  the  length  of  the  separated  region,  with  the  move¬ 
ment  of  the  separation  point  on  to  the  curved  segment.  The  distributions  to  the 
configurations  with  ramp  angles  of  24*  {shown  in  Fig.  9b)  and  30*  were  prac¬ 
tically  identical  to  those  obtained  on  the  flat  plate -wedge  models  with  the  same 
ramp  angles  (see  Fig.  8c). 

The  first  measurable  effect  of  introducing  surface  curvature  was  evi¬ 
dent  when  the  flat  plate- 12"  radius  section- 18”  ramp  was  tested.  The  distri¬ 
butions  shown  in  Fig.  10a  reveal  a  rather  important  feature.  The  skin  friction 
distribution  shows  that  a  small  region  of  separated  flow  exists  on  the  curved 
ramp.  The  pressure  and  heat  transfer  measurements  do  not  give  clear  evidence 
of  this.  Although  the  pressure  distribution  displays  some  evidence  of  an  inflex¬ 
ion  point  in  the  separated  region,  neither  it  nor  the  local  heat  transfer 
distribution  could  be  cited  to  point  out  the  existence  of  a  separated  region. 

Here  and  in  the  measurements  of  separated  regions  in  entropy  layers  gener¬ 
ated  by  strong  bluntness,  reported  in  Ref.  29,  skin  friction  measurements 
are  needed  to  determine  accurately  when  separation  first  occurs.  Adding  a 
6°  segment  to  turn  the  flow  through  24*  causes  a  large  separated  region  with 
well  defined  separation,  plateau  and  reattachinent  regions  in  the  pressure 
distribution  (see  Fig .  10b). 
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The  distributions  obtained  over  two  of  the  configurations  with  the  20" 
radius  of  cylindrical  section  are  shown  in  Figs.  11a  and  lib.  For  this  sur¬ 
face  curvature  the  wedge  angle  for  incipient  separation  has  been  increased  to 
approximately  20*;  however,  when  the  wedge  angle  is  increased,  the  size  of 
the  separated  region  increases  very  rapidly.  The  distribution  for  the  24* 
ramp  angle,  shown  in  Fig.  11a,  shows  the  separated  region,  the  leading  edge 
of  which  has  fed  forward,  covering  the  curved  surface.  Once  the  separated 
region  moves  upstream  of  the  flat  plate -curved  surface  junction,  then  the 
subsequent  development  and  structure  for  increasing  wedge  angle  is  similar 
to  that  observed  on  a  flat  plate-wedge  model  (see  Fig.  lib). 

Finally,  two  of  the  distributions  on  the  12"  and  20"  radius  compression 
surfaces  are  shown  in  Figs.  12a  and  12b  for  the  highest  unit  Reynolds  number 
conditions.  Increasing  the  unit  Reynolds  number  increases  the  ratio  of  the 
radius  of  curvature  to  the  local  boundary  layer  thickness,  which,  if  curvature 
effects  are  important,  should  cause  a  decrease  in  the  length  of  the  separated 
region;  while  we  know  from  the  measurements  on  flat  plate-wedge  surfaces 
that  increasing  the  unit  Reynolds  number  causes  an  increase  in  the  size  of  the 
separated  region  or  brings  an  attached  region  closer  to  separation.  In  each 
of  the  cases  shown  here,  increasing  the  unit  Reynolds  number  moved  the  separa¬ 
tion  point  further  forward.  The  curvature  of  the  surface  appeared  to  be  of 
little  importance  once  a  separated  region  was  established. 

3.5  PITOT  AND  CONE  PRESSURE  MEASUREMENTS 

Measurements  of  the  pitot  and  cone  pressure  were  made  across  the 
boundary  layer  in  the  separation  and  reattachment  regions  on  the  models  in 
order  to  obtain  an  estimate  of  the  pressure  gradients  normal  to  the  surface 
in  these  regions.  A  photograph  of  the  rake  containing  the  conical  probes  is 
shown  in  Fig.  7.  The  probes,  which  are  1/8-in  diameter,  have  the  piezoelectric 
sensing  element  and  the  F.E.T.  impedance  transform  circuit  mounted  on  the 
cylindrical  section  close  to  the  orifices.  The  conical  tip  and  the  pitot  head 
are  interchangeable  so  that  the  same  sensing  element  is  used  to  record  both 
cone  and  pito.  pressures.  The  Mach  number  at  the  tip  of  the  probe  can  be 
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obtained  from  the  ratio  of  the  pitot  and  cone  pressures;  then  from  the  Mach 
number  and  cone  pressure,  the  local  static  pressure  can  be  estimated. 
Although  this  technique  is  accurate  at  low  supersonic  Mach  numbers,  because 
the  ratio  of  cone  to  total  pressure  is  less  sensitive  to  Mach  number  in  hyper¬ 
sonic  flow,  the  accuracy  is  estimated  to  be  no  better  than  20%  over  the  range 
of  conditions  at  which  our  measurements  were  made.  Fig.  13  shows  schlieren 
photographs  and  the  derived  static  pressure  distributions  for  measurements 
made  in  the  separation  and  reattachment  regions.  The  measurements  made 
just  downstream  of  separation  do  not  show  the  presence  of  a  strong  pressure 
gradient  normal  to  the  surface.  However,  the  measurements  in  and  above  the 
reattaching  boundary  layer  indicate  that  the  normal  pressure  gradient  in  this 
region  is  comparable  with  that  in  the  streamwise  direction.  Although  the 
exact  height  of  the  boundary  layer  cannot  be  determined,  using  an  estimate 
from  the  schlieren  photographs,  the  total  static  pressure  ratio  across  the 
boundary  layer  is  found  to  be  approximately  1.9.  This  figure  is  of  the  same 
magnitude  as  predicted  from  theory;  however,  it  must  be  emphasized  that  the 
probe  measurements  are  intended  to  serve  only  as  a  qualitative  estimate. 

3 . 6  INCIPIENT  SEPARATION 

In  this  study  we  did  not  find  the  strong  dependence  of  incipient  separa- 
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tion  on  surface  curvature  that  was  observed  in  the  earlier  studies  of  Kuehn, 
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and  Sterrett  and  Emery.  A  correlation  of  the  incipient  conditions  obtained 
in  the  present  study  is  shown  in  Fig.  14.  Here  we  have  plotted  the  product  of 
the  free  stream  Mach  number  and  the  angle  to  promote  incipient  separation, 

Moo  ®map  »  versus  ,  the  product  of  the  viscous  interaction 

parameter  and  the  ratio  of  the  Reynolds  number  based  on  the  radius  of  curva¬ 
ture  of  the  circular  insert  and  the  local  boundary  layer  thickness.  Earlier 
studies  had  shown  that  for  sharp  flat  plate-wedge  configurations  M Q,nc,p  **-  • 

and  for  a  fixed  radius  of  curvature  the  configurations  tested  have  also  exhibited 
this  trend.  The  ratios  of  radius  of  surface  curvature  to  boundary  layer  thick¬ 
ness  varied  from  0  to  22,  similar  to  those  obtained  by  Kuehn,  et  al;  however, 
only  at  the  largest  values  of  this  ratio  was  a  significant  increase  found  in  the 
wedge  angle  to  promote  incipient  separation.  It  must  be  noted  that  conditions 
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encountered  in  the  present  study  differ  greatly  from  those  obtained  by  Kuehn, 
and  Sterret  and  Emery.  Here  our  experiments  were  conducted  at  Mach  num¬ 
bers  which  were  2  to  5  times  larger  which,  from  the  earlier  studies,  should 
amplify  the  curvature  effect.  However,  the  boundary  layers  in  the  present 
study  were  completely  laminar  as  compared  with  the  transitional  or  fully 
turbulent  conditions  encountered  earlier,  and  the  measurements  were  made 
on  a  highly  cooled  wall,  whereas  Kuehn's  an*  Sterrett  and  Emery's  experi¬ 
ments  were  conducted  over  adiabatic  surfaces. 


We  believe  that  relative  insensitivity  of  the  separation  condition  to 
surface  curvature  results  from  the  intrinsic  structure  of  the  laminar  boundary 
lave  r  on  highly  cooled  high  Mach  number  flows.  In  these  flows  most  of  the 
momentum  contained  in  the  boundary  layer  is  concentrated  at  its  outer  edge; 
thus,  in  order  to  produce  a  change  in  its  development,  a  slope  change  at  the 
surface  must  be  propagated  almost  completely  across  the  boundary  layer.  If 
the  distance  travelled  downstream  by  this  disturbance  (  *  Mf  S  )  is  small 
compared  with  the  projected  distance  along  the  edge  of  the  boundary  layer  to 
the  compression  surface,  then  the  initial  curvature  of  surface  is  important. 
However,  if,  as  in  the  cases  studied  here,  the  boundary  layer  begins  to  turn 
parallel  to  the  curved  surface  only  when  it  is  in  close  proximity  to  this  surface, 
then  it  is  the  surface  curvature  at  this  point  which  is  of  importance.  The  angle 
which  the  edge  of  the  viscous  layer  makes  with  the  local  angle  of  the  ramp 
just  before  appreciable  turning  takes  place  is  of  principal  importance.  If 


this  angle  is  less  than  or  equal  to  the  wedge  angle  for  incipient  separation, 
in  ip  *  4^en  we  w°uld  not  expect  separation  to  occur.  Following  this 
reasoning,  a  significant  change  could  be  expected  in  the  incipient  separation 
condition  by  introducing  surface  curvature  only  if - -incip  >  .  Since 

_  t/g  ^ 

M0W  oC  t  ,  this  criteria  implies  that  at  higher  Mach  numbers,  for  a 

map  ll 

fixed  length  of  flat  plate,  an  even  larger  radius  is  required  before  surface 


curvature  effects  are  of  importance. 
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4. 


CONCLUDING  REMARKS 


In  this  paper  a  theoretical  analysis  has  been  developed  to  describe  the 
non- similar  development  o£  attached  and  separated  laminar  boundary  layers 
in  regions  of  shock,  wave -boundary  layer  interaction  over  highly  cooled  curved 
compression  surfaces.  No  assumption  was  made  that  the  pressure  is  constant 
across  the  boundary  layer;  thus,  by  obtaining  a  simultaneous  solution  to  the 
integral  forms  of  the  conservation  equations  including  the  normal  momentum 
equation,  we  were  able  to  calculate  the  pressure  difference  across  the  bound¬ 
ary  layer  and  the  effect  of  surface  curvature.  The  initial  response  of  a 
boundary  layer  to  a  streamwise  pressure  gradient  depends  strongly  on  the 

normal  pressure  gradient;  small  values  of  the  parameter  of  — *  ~ ..ft*  changed 

A* 

the  boundary  layer  from  supercritical  to  subcritical.  The  analysis  was  com¬ 
pared  with  measurements  reported  in  this  and  in  earlier  papers  and  was 
found  to  be  in  good  agreement.  These  calculations  indicated  that  pressure 
ratios (^~)  of  the  order  of  1 . 6  are  generated  in  the  separation  and  reattach¬ 
ment  regions  of  the  flow. 

In  the  experimental  study,  skin  friction,  heat  transfer  and  pressure 
measurements  were  made  on  a  series  of  flat  plate -cylindrical  arc-wedge 
compression  surfaces  to  examine  the  effects  of  radius  of  curvature  on  the 
occurrence  and  structure  of  separated  regions.  Pitot  and  cone  pressure 
traverses  were  also  made  to  obtain  an  estimate  of  the  pressure  gradients 
normal  to  the  surface  in  the  separation  and  reattachment  regions.  In  these 
flows,  the  first  indication  that  flow  separation  was  to  occur  was  provided  by 
the  skin  friction  measurements  on  the  compression  surface;  the  skin  friction 
was  first  observed  to  be  negative  on  the  surface  downstream  of  the  flat  plate. 
On  the  compression  surfaces  with  large  radii  of  curvature,  the  incipient 
separa??un  condition  could  not  be  detected  with  any  accuracy  by  observing  the 
first  occurrence  of  an  inflexion  point  in  the  pressure  distribution  or  a  change 
in  the  form  of  the  heat  transfer  distribution.  The  static  pressure  distribu¬ 
tions,  evaluated  from  the  pitot  and  cone  pressure  measurements  across  the 
boundary  layers  in  the  separation  and  reattachment  regions,  indicated  that 
the  normal  pressure  gradients  were  of  the  same  magnitude  as  those  in  the 


27 


streamwise  direction.  Although  the  experiments  were  conducted  for  a  range 
of  ratios  of  S0/R  which  were  similar  to  those  used  by  Kuehn  and  Sterrett 
and  Emery,  we  found  that  incipient  separation  was  far  less  sensitive  to  sur¬ 
face  curvature  than  indicated  in  these  earlier  studies.  This  we  believe  results 
from  the  structure  of  the  laminar  boundary  layer  in  hypersonic  flow  over 
highly  cooled  walls  where  the  streamwise  distances  for  a  surface  disturbance 
to  propagate  to  the  edge  of  the  boundary  layer  (where  most  of  the  momentum 
is  concentrated)  are  greater  than  the  length  of  the  curved  surface. 
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APPENDIX  1 


The  boundary  layer  equation  can  be  written  in  the  form 


>  **  °9'  5‘  */dx  lo3*  ax  CV  «rx  = 


Continuity  Equation  =  f  ) 
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X-  Momentum  Equation  (^  -  2  ) 
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x- Moment  of  Momentum  Equation  (j*  J  ) 


tf- Momentum  Equation  (  j  f  ) 
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Energy  Equation  (  j -  5  ) 
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evaluation  of  the  integral  functions 


Cohen  and  Reshotko  Profiles 


The  integral  functions  H ,  J  ,  W  ,  Q  ,  p  and  R  were  evaluated 
from  the  family  of  velocity  profiles  derived  by  Cohen  and  Reshotko  and 
expressed  in  the  form 


H  -  ca  +  c ,  a.  y-  c2a* 


c  a 


Similarly!  integral  functions  which  depended  only  upon  the  enthalpy  profiles 
were  fitted  as  a  power  series  in  b  .  The  functions  1~  ,  t  ,  which  are 
functions  of  both  a,  and  6  ,  we  have  written  in  the  form 


1  sw 


cLY 


J>k) 


I  (b) 

G«a) 


T  r 


Pr  *  - 


<4 


cLY 


K(*.  b  ) 

G  (a.) 


where  K  has  been  expressed  as  a  function  of  both  cl  and  fa  in  the  form 
Km  -  c.  *  c.  b  *  c  A  ’  +  c.  6  3 . 

4L~  r  1r  2r  Jr 

The  numerical  values  for  the  coefficients  in  these  power  series  are  presented 
in  Ref.  28. 


Polynomial  Profiles 

The  integral  functions  P f  Pz  .  .  .  Pl9  were  evaluated  by  exact 
integration  of  the  quartic  profiles  where  the  Cohen  and  Reshotko  and  poly¬ 
nomial  profiles  are  connected  by  the  relationships. 

For  attached  flow  A  =  a.  and  for  separated  flow  A-  (f-a.). 

For  attached  flow 

8  =  -12.045302  +  b (177. 25491  +  b (-942.  39173  +  b  (2530. 4358 
+  6(-3382.  3794  +  b  1782.0275)))) 

For  separated  flow 

8  =  0.024491381  +  b(10.  530295  +  b(-27. 018912  +  b30.  097469)) 

The  values  of  Pf  •  •  •  P,5  are  given  below . 
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Figure  3  A  COMPARISON  BETWEEN  COHEN  AND  RESHOTKO  AND  POLYNOMINAL  PROFILES 


Figure 4  DISTRIBUTION  OF  PROFILE  PARAMETERS 
(M^  -  10,  ReL  -  8.1  x  105, 8^  «  12.5°) 
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Figure  5  A  COMPARISON  BETWEEN  THEORY  AND  MEASUREMENTS 
OF  PRESSURE  AND  HEAT  TRANSFER 

(Mqo=  10,  ReL  =  8.1  x  105,  0W  =  12.5°) 
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Figure  6b  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENTAL 
PRESSURE  MEASUREMENTS 

(Mqo=  14.5,  Re/FT=  1.8  x  105,  8W  =  18°) 
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Figure  6c  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENTAL 
HEAT  TRANSFER  MEASUREMENT 

(M  =  14.5,  Re/FT  =  1.8  x  105,  0W  =  18°) 
oo  w 
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c.  WEDGE  ANGLE  24° 

Figure  8  SKIN  FRICTION,  HEAT  TRANSFER,  AND  PRESSURE  DISTRIBUTIONS 
ON  THE  FLAT  PLATE  -  WEDGE  MODELS 

(M„  -  14.1,  Re/FT-  7.2  x  104) 
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b.  WEDGE  ANGLE  24° 

Figure  9  SKIN  FRICTION,  HEAT  TRANSFER  AND  PRESSURE  DISTRIBUTIONS 
ON  THE  FLAT  PLATE  -  6"R  CYLINDRICAL  ARC-WEDGE  MODELS 

(M.  -  14.1,  Re/FT  =  7.2  x  104) 
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Figure  10  SKIN  FRICTION,  HEAT  TRANSFER  AND  PRESSURE  DISTRIBUTIONS 
ON  THE  FLAT  PLATE  -  12"R  CYLINDRICAL  ARC-WEDGE  MODELS 

<M„  -  14.0,  R«/FT  -  7.0  x  104) 


50 


b.  WEDGE  ANGLE  30° 


Figure  1 1  SKIN  FRICTION.  HEAT  TRANSFER  AND  PRESSURE  DISTRIBUTIONS 
ON  THE  FLAT  PLATE  -  20"R  CYLINDRICAL  ARC-WEDGE  MODELS 
(M^  *  14.0,  Re/FT  -  7.2  x  104) 
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Figure  12  SKIN  FRICTION.  HEAT  TRANSFER  AND  PRESSURE  DISTRIBUTIONS 
ON  THE  FLAT  PLATE  -  CYLINDRICAL  ARC-WEDGE  MODELS 
(M_  -  14.5,  Re/FT  -  1.8  x  108) 


Figure  13  STATIC  PRESSURE  MEASUREMENTS  ACROSS  THE  BOUNDARY  LAYER 

(Mm  =  14.5,  Re/FT  =  1.8  x  105,  0  =  18°) 
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